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The ground-state polarization of PbTiC>3 thin films is studied using a microscopic effective 
Hamiltonian with parameters obtained from first-principles calculations. Under short-circuit elec- 
trical boundary conditions, (001) films with thickness as low as three unit cells are found to have a 
perpendicularly polarized ferroelectric ground state with significant enhancement of the polarization 
5— l ' at the surface. 

' Ferroelectricity is a collective phenomenon and, as such, it is expected to be strongly influenced by surfaces and 
finite size effects . It has been generally thought that decreasing the thickness of thin films of perovskite ABO3 
■ compounds suppresses ferroelectricity and eliminates it altogether at a nonzero critical thickness. In PbTi03, this 
critical size has been estimated as 7-10 nm However, recently synthesized high quality films of PbZr^Tii^^Oa 

1 ^ '■ have been observed to exhibit a stable perpendicularly polarized ferroelectricity down to and possibly even below 
a thickness of 40 A ||. This suggests that the earlier observations could stem from variations in electrical and 
1 ' mechanical boundary conditions or from extrinsic effects associated with defects and impurities. 

In this Letter, we address the question of the intrinsic ferroelectricity of PbTiC>3 thin films using a first-principles 

S effective Hamiltonian CHeff) based on that previously developed for studies of the finite-temperature ferroelectric 
transition in the bulk [|7|. This approach combines the advantages of phenomenological models, which are physically 
transparent, with those of first-principles simulations, which yield highly accurate material-specific microscopic infor- 
mation. We find that PbTiC>3 films exhibit a ferroelectric instability which is even greater than that of the bulk, and 
identify the reason for this within the model. Further, we analyze the spatial variation of the polarization within the 
film, showing that polarization is significantly enhanced over the bulk value in a single surface layer. 

We consider symmetrically TiC>2 terminated (001) slabs of PbTiC>3, under stress-free and short-circuit boundary 
conditions. The latter are imposed by the presence of perfectly conducting sheets positioned at a distance a /4 above 
and below the top and bottom surfaces of the slab, where ao is the experimental lattice constant of bulk PbTiC>3 
(a =3.97 A). 

As for the bulk 0, for the thin films, Tteff is constructed by projecting the full interionic Hamiltonian into a 
^-j- . subspace which contains the relevant degrees of freedom. These are the localized atomic displacement patterns (&) 
' centered around each unit cell i, determined by the lattice Wannier function method |^2| , and the homogeneous strain 
C*") . tensor e a p. The form of the resulting model is 
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^~^eff — (^-dipolar H~ ^sr) H~ ^-anharm 



where the harmonic part is separated into a long-range dipolar term, Hdipoiar, and short-range corrections, 7i sr . 
T~Lanharm represents the anharmonic couplings and distrain combines the elastic energy and the lowest order coupling 
of £i to e a p. The local dipole per unit cell is pi = Z*ea^i where Z* is the Born effective charge associated with £j. 

The present construction, based on projection from a thin-film interionic Hamiltonian, is substantially refined over 
that presented previously jg], which was a simple slab truncation of the bulk H. e ff- However, in order to avoid 
computationally intensive first-principles slab calculations |p[-pl|, we use an approximate interionic Hamiltonian, 
obtained by transfer of microscopic information available for the bulk to the thin film geometry. The associated 
modifications include termination of the interatomic short-range force constants at the surface |13|] , change in the 
effective dipole-dipole interaction resulting from the perfectly conducting plates which implement the short-circuit 
boundary conditions, and corrections to preserve global translational symmetry and charge neutrality fll4| . The 
projection is then performed using the same basis functions as for the bulk, with minor modifications for the unit cells 
at the surface. The anharmonic and strain terms are, for simplicity, assumed to be the same as in the bulk. Further 
details will be given in Ref. Jji| . 

The minimum of Ti e f f for in-plane and perpendicularly polarized films is reported for different thicknesses in Table Q. 
The polarization is assumed to be uniform within each layer. The mechanical boundary conditions are taken as stress 
free, so that the strain relaxes to its optimum value. At each thickness, perpendicular polarization is seen to be most 
favorable, with stabilization relative to the paraelectric state increasing with decreasing thickness. The corresponding 
layer-by-layer polarization profiles are shown in Fig. [l]. While in the interior the polarization approaches the bulk 
value even for very thin films, it is significantly enhanced at the surface. In fact, unpolarizcd films are unstable 
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against polarization of only the two surface layers. The enhancement of the average polarization couples to the strain, 
resulting in a c/a ratio that increases slightly with decreasing thickness. 

To understand the ferroelectric instability and surface polarization enhancement in perpendicularly polarized films, 
we consider Ti. e f t specialized to the case of perpendicular polarization, constant in each layer. The computed param- 
eters are both layer and thickness dependent, approaching bulk values in the interior of the film. As illustrated in 
Table I and in Figure [l] this behavior can be closely reproduced by assigning bulk values to all parameters except 
the onsite interaction in the surface layer and the quadratic interaction between the surface and its neighboring layer, 
and by approximating the dipolar interaction matrix as described below. This corresponds to the following simplified 
model with parameters reported in Table [n| The short-range and anharmonic part take the form: 

JV N-l AT-2 

H sr + n anharm = ]T(a, C + bC + c£ + d C) + E a + E b Z» 

n— 1 71 — 1 71— 1 

+AA t (C 2 + &) + Aa 6 

where Ai, B, C, D, a and b are directly related to bulk tieff parameters, while AAi and Aa arise from changes at the 
surfaces and are, to a good approximation, independent of the thickness. 

The dipolar contribution, computed using the Ewald summation technique including the effect of the metallic plates, 
is accurately represented (to within 0.1%) by the expression 
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Udipolar - ~ [ , AT C °T/n\ (52 Z n^i) 2 + U se if V] Z* 2 £ (1) 
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where the optical dielectric constant, Eoo, is taken equal to 8.24 ||7|,[L3| an d Z* is +10.218 at the surface and +10.056 
elsewhere. 

Finally, the terms which represent the coupling of the polarization to homogeneous strain are derived directly from 
the corresponding bulk terms: 

N N 

T-tstrain = zl ^ + NCl2 X/ e <=*<± e PP + (zJ ^n)(p0 2_j Gaa + 91^33.) ■ 

ol a<{3 n—1 a 

Within this simplified model, the origin of the enhancement of the ferroelectric instability in thin films over the bulk 
can be readily identified. The model can be obtained from a reference model which has the uniform bulk ferroelectric 
ground state by modifying three terms as follows: (i) changing the dipolar contribution by moving the perfectly 
conducting plates, at the level of the model, outward by ao/4 from where the array of dipoles combined with their 
images forms an infinite simple cubic lattice; (ii) suppressing the periodic boundary conditions on the short-range 
terms, that is, eliminating the coupling of layers 1 and 2 to layers N and N+l , and (iii) including the short-range 
surface corrections AA[ and Aa. We can introduce these modifications of the bulk model separately. 

By examining first the effects of term (i), we can qualitatively understand how the electrical boundary condition 
plays a crucial role in determining the thin film ground state [fl6|| . The change in the dipolar contribution from the 
bulk to the thin film significantly suppresses the perpendicular polarization. This is not surprising as motion of the 
plates outward from the slab decreases the compensation of the depolarization field. In contrast, starting from the 
bulk and truncating the antiferroelectric short-range couplings at the surface (term (ii)) leads to a large energy gain, 
a large enhancement of the surface polarization, and a somewhat lesser enhancement of the polarization of the interior 
layers which is close to constant starting with the second layer. Examining the effects of the short-range surface terms 
A Ai and Aa (term (iii)), we find similar behavior, except that the polarization in the second layer, while still greater 
than its bulk value, is suppressed relative to that of the interior layers. 

In the full thin film model, the combined effects of terms (ii) and (iii) compete with the suppression of ferroelectricity 
by the dipolar contribution of term (i) . The surface enhancement and characteristic shape of the polarization profile 
is preserved by the inclusion of term (i) , while the magnitude of the overall enhancement is lowered, so that the values 
of polarization in the interior approach the bulk value. 

When comparing to observations on ultrathin PbTiC>3 films, several factors must be taken into account. First, 
substrates will generally induce epitaxial stress, which can significantly change the ferroelectric state of the film p7|. 
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The screening of the depolarization field by surface charges can be expected to be less effective than for the ideal case 
of perfectly conducting plates considered here, which would act to reduce the polarization enhancement Jl8| . Finally, 
the finite conductivity of PbTi03 could lead to the creation of a depletion layer and related effects jl9| not included 
in this model. Nevertheless, the plausibility of surface polarization enhancement warrants further investigation. 

In conclusion, we have constructed a microscopic H e f / for the study of PbTiC>3 thin films under stress- free and 
short-circuit boundary conditions. This model shows that (001) PbTiG<3 films as thin as three unit cells exhibit a 
perpendicularly polarized ferroelectric ground state, with significant enhancement of the polarization at the surface. 
The ferroelectric instability is consistent with recent observations, suggesting further theoretical and experimental 
investigation. 
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FIG. 1. Layer-by-layer profiles of the polarization for films of thickness N=3 (squares), N=5 (solid circles) and N=7 
(diamonds). The polarization profile for the simplified model for N=5 is shown by the open circles and dashed line. For the 
bulk, the computed polarization in the tetragonal phase is equal to 0.83 C/m 2 . 



TABLE I. Energy per unit cell in eV (in reference to the paraelectric state) and macroscopic strain as a function of thickness 
(N) for perpendicularly (_L) and in-plane (||) polarized films. The last column refers to the simplified model (see text). For 
comparison, in the tetragonal bulk phase: E = —0.086 eV, e xx = e yy = —0.015 and e zz = 0.064. 
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3 


-0.110 


-0.017 


0.071 


-0.320 


-0.023 


0.098 


-0.309 


5 


-0.076 


-0.014 


0.061 


-0.217 


-0.020 


0.084 


-0.210 


7 


-0.072 


-0.014 


0.060 


-0.163 


-0.018 


0.076 


-0.172 



TABLE II. Parameters of the simplified Tieff when energies are in eV, lengths in A and £ n is normalized to ao- 

A t -66.2664 a +118.046 u conv -14.3997 Cn +117.9 

B +6.906 x 10 3 b +14.718 u self +32.7745 C12 +51.5 

C -1.658 x 10 5 AAi -46.866 u surf -0.0472 g -107.7 

D +9.630 x 10 6 Aa +26.598 u nn -2.3788 gi -790.3 
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